Accelerating Cold Dark Matter Cosmology (Qa = 0) 
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Abstract 

A new kind of accelerating flat model with no dark energy that is fully dominated by cold dark matter 
(CDM) is investigated. The number of CDM particles is not conserved and the present accelerating 
stage is a consequence of the negative pressure describing the irreversible process of gravitational 
particle creation. A related work involving accelerating CDM cosmology has been discussed before 
the SNe observations [Lima, Abramo & Germano, Phys. Rev. D53, 4287 (1996)]. However, in order 
to have a transition from a decelerating to an accelerating regime at low redshifts, the matter creation 
rate proposed here includes a constant term of the order of the Hubble parameter. In this case, Hq 
does not need to be small in order to solve the age problem and the transition happens even if the 
matter creation is negligible during the radiation and part of the matter dominated phase. Therefore, 
instead of the vacuum dominance at redshifts of the order of a few, the present accelerating stage in 
this sort of Einstein-de Sitter CDM cosmology is a consequence of the gravitational particle creation 
process. As an extra bonus, in the present scenario does not exist the coincidence problem that 
plagues models with dominance of dark energy. The model is able to harmonize a CDM picture 
with the present age of the universe, the latest measurements of the Hubble parameter and the 
Supernovae observations. 

PACS numbers: 98.80.-k, 95.35.+d,95.30.Tg 



I. INTRODUCTION 

A large amount of data relevant to cosmology (involv- 
ing Supernovae type la and cosmic background radia- 
tion probes) have provided strong evidence that the ob- 
served universe is well described by an accelerating, flat 
Friedmann- Robertson- Walker (FRW) model P, H, 11, |1|- 
However, the substance or mechanism behind the current 
cosmic acceleration remains unknown and constitutes a 
challenging problem of modern cosmology. 

In relativistic cosmology, an accelerating regime is ob- 
tained by assuming the existence of a dark energy com- 
ponent (in addition to cold dark matter), an exotic fluid 
endowed with negative pressure in order to violate the 
strong energy condition Q . The simplest theoretical rep- 
resentation of dark energy is by means of a cosmological 
constant A, which acts on the Einstein field equations 
(EFE) as an isotropic and homogeneous source with con- 
stant equation of state (EoS) w = p/ p = — 1. 

All observational data available so far seems to be in 
good agreement with the cosmic concordance model, i.e., 
a vacuum energy plus cold dark matter (ACDM) sce- 
nario. Nevertheless, ACDM models are plagued with 
several problems. For instance, it is very difficult to 
reconcile the small value required by observations (~ 
10^^''erg/cm'^) with estimates from quantum field the- 
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ories ranging from 50-120 orders of magnitude larger [6(]. 
Such problem has inspired many authors to propose alter- 
native candidates in the literature [3, [1, |^, [l3, [HI , among 
them: (i) a relic scalar field slowly rolling down its poten- 
tial, (ii) a A(i)-term or a decaying vacuum energy density, 
(iii) the "X- matter" , an extra component characterized 
by equation of state Px — i^Px , where oj may be constant 
or a redshift dependent function, (iv) a Chaplygin-type 
gas whose equation of state is p = —A/p°', where A and 
a are positive parameters. More recently, some attention 
has also been paid to a possible interaction between the 
dark sector components [l^l . 

The space parameter of such models are usually highly 
degenerated and some of them contain the ACDM sce- 
nario as a particular case. In point of fact, the plethora of 
possible candidates does not help to identify the nature 
of this mysterious component since there is no compelling 
direct evidence yet for dark energy (or its dynamical ef- 
fects). In other words, the evidence supporting its exis- 
tence is not strong enough to be considered established 
beyond doubt (see [l^ for a critical discussion). 

Roughly speaking, a realistic cosmological scenario 
should be in agreement with at least four well estab- 
lished observational results, namely: (i) the existence of 
a dark non-baryonic component as required by the dy- 
namics of galaxies and clusters, the matter power spec- 
trum and other independent probes like the temperature 
anisotropies of the cosmic microwave background (CMB) 
from last scattering surface (ii) the late time cosmic ac- 
celeration, (iii) the (nearly) flatness of the Universe, and, 
finally, (iv) a Hubble parameter i/o ~ 72 km/s.Mpc with 
the Universe being older than 12 Gyrs in order to accom- 



modate the oldest observed structures (globular clusters). 
When confronted with this simple requirements, we see 
that the CDM or Einstein-de Sitter cosmology is in clear 
contradiction with results (ii) and (iv). Therefore, if one 
assumes that the dark energy does not exist, the first 
task is to explain how a flat CDM dominated Universe 
can accelerate at late times because, potentially, acceler- 
ating cosmologies solve the age problem. 

In this concern, we recall that the presence of a nega- 
tive pressure is the key ingredient required to accelerate 
the expansion. This kind of stress occurs naturally in 
many different contexts when the physical systems de- 
part from a thermodynamic equilibrium states |14l |. In 
general, such states are connected with phase transitions 
(for example, in an overheated van der Waals liquid), and 
for some systems the existence of negative pressure seems 
to be inevitable llSll . In this connection, as first pointed 
out by Zeldovich [1^ , the process of cosmological particle 
creation at the expenses of the gravitational field can phe- 
nomenologically be described by a negative pressure and 
the associated entropy production. In principle, such an 
approach is completely different from the one developed 
by Hoyle and Narlikar [l3| adding extra terms to the 
Einstcin-Hilbert action describing the so-called C-field. 
In the latter case, the creation phenomenon is explained 
trough a process of interchange of energy and momen- 
tum between matter itself and the C-field as happens, 
for instance, in vacuum decaying cosmologies |8|. 

The gravitational matter creation processes was in- 
vestigated frorn^a microscopic viewpoint by Parker and 
collaborators |18| by considering the Bogoliubov mode- 
mixing technique in the context of quantum field theory 
in curved space-time [131 . Despite being rigorous and 
well-motivated, those models were never fully realized, 
probably due to the lack of a well-defined prescription of 
how matter creation is to be incorporated in the classical 
EFE. 

The consequences of gravitational matter creation 
have also been macroscopically investigated mainly as 
a byproduct of bulk viscosity processes near the Planck 
era as well as during the reheating of the infiationary sce- 
narios [2y, [2l|. However, the first self-consistent macro- 
scopic formulation of the matter creation process was put 
forward by Prigogine and coworkers f22] and somewhat 
clarified by Calvao, Lima and Waga [23] through a mani- 
festly covariant formulation. It was also shown that mat- 
ter creation, at the expenses of the gravitational field, can 
effectively be discussed in the realm of the relativistic 
nonequilibrium thermodynamics. Later on, it was also 
demonstrated that the matter creation is an irreversible 
process completely different from bulk viscosity mecha- 
nism [2^ (see also [28| for a more complete discussion). 
Several interesting features of cosmologies with creation 
of matter and radiation have been investi gate d by many 
authors [H, [H, [13, IH, ^, [U (see also [Hi for recent 
studies on this subject). 

In comparison to the standard equilibrium equations, 
the irreversible creation process is described by two new 



ingredients: a balance equation for the particle number 
density and a negative pressure term in the stress tensor. 
Such quantities are related to each other in a very definite 
way by the second law of thermodynamics [22, [23l ■ The 
leitmotiv of this approach is that the matter creation pro- 
cess, at the expense of the gravitational field, can happen 
only as an irreversible process constrained by the usual 
requirements of non-equilibrium thermodynamics. 

In this context, we are proposing here a new flat cosmo- 
logical scenario where the cosmic acceleration is powered 
uniquely by the creation of cold dark matter particles. 
As we shall see, the model is consistent with the super- 
novae type la data, and a transition redshift of the order 
of a few is also obtained. In this extended CDM model, 
the Hubble parameter does not need to be small in order 
to solve the age problem and the transition happens even 
if the matter creation is negligible during the radiation 
and considerable part of the matter dominated phase. 
Moreover, the so-called coincidence problem of dark en- 
ergy models is replaced here by a gravitational particle 
creation process at low redshifts. 



II. COSMOLOGY AND MATTER CREATION 

For the sake of generality, let us start with the homo- 
geneous and isotropic FRW line element 



ds^ =dt^ -R^{t){ 



dr"^ 



1-kr 



r^d0^ + r^sin^9d(j}^), (1) 



where R is the scale factor and fc = 0, ±1 is the curvature 
parameter. Throughout we use units such that c — 1. 

In that background, the nontrivial EFE for a fluid 
endowed with matter creation and the balance equa- 
tion for the particle number density can be written as 
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where an overdot means time derivative and p, p, Pc, n 
and V' are the energy density, thermostatic pressure, cre- 
ation pressure, particle number density and matter cre- 
ation rate, respectively. The quantity F with dimension 
of {time)~^ is the creation rate of the process. The cre- 
ation pressure Pc is defined in terms of the creation rate 
and other physical quantities. In the case of adiabatic 
matter creation, it is given by [22, [23, [2J, [2^, [2^ (see 
also Appendix A for a simplified deduction) 
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where H = R/R is the Hubble parameter. 

As one may check, by combining the EFE with usual 
equation of state, p = cup, the equation governing the 
evolution of the scale function is readily obtained: 
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The above expression shows how the matter creation 
rate, F, modifies the evolution of the scale factor as com- 
pared to the case with no creation. Conversely, the cos- 
mological dynamics with irreversible matter creation will 
be defined once the matter creation rate is given. As 
should be expected, by taking F = it reduces to the 
FRW differential equation governing the evolution of a 
perfect simple fluid |33| . 




FIG. 1: The scale function as a function of time. Like the 
FRW dust filled model (solid black line), the evolution starts 
from an initial singularity. Note that the influence of 7 is 
negligible in the early universe (see Eq. (|12|l for H ^ -ffo)- 



III. FLAT CDM MODEL WITH MATTER 
CREATION AND THE AGE OF THE UNIVERSE 

In what follows we focus our attention on the flat cold 
dark matter model {k ^ lu = 0) with the previous equa- 
tion reducing to: 
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F 
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R^ = 0, 



or, equivalently. 



H+-H-' 1--— =0. 



2 V 3iJ, 
On the other hand, Eq. (j4|) can be rewritten as 
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which means that the creation process can effectively be 
quantified by the dimcnsionless ratio (see also Eq. ([5])) 



A{t) 



F 
3^' 



(10) 



which in general is a function of time. If F ^ 3H, that 
is, A <C 1, the creation process is negligible leading to 
n oc i?^"^ and H = 2/3t, as should be expected for an 
Einstein-de Sitter model. The opposite regime (F 3> 3H) 
defines an extreme theoretical situation, where the cre- 
ation process is a phenomenon so powerful that the dilu- 
tion due to expansion is more than compensated. Prob- 
ably, this kind of behavior may happen only in the very 
early universe as happens, for instance, during the re- 
heating stage of inflation. An intermediary (and phys- 
ically more reasonable situation) occurs if this ratio is 



smaller or of the order of unity (F < 'iH). In particu- 
lar, if F = 3H the dilution due to expansion is exactly 
compensated and the number density remains constant. 
From now on we consider that A(i) < 1. 

In a series of papers [2y, [23], we have investigated 
some properties of adiabatic matter creation models with 
F — 3f3H, where /3 is a constant parameter contained on 
the interval [0,1] (A — (3). However, that kind of models 
are always accelerating for /? > 1/3 or decelerating for 
(3 < 1/3, that is, there is not a transition redshift from a 
decelerating to an accelerating regime as required by the 
supernovae type la observations (see Figure 3a). In order 
to cure such a difficulty we add a constant term in this 
expression, that is, we consider the following matter cre- 
ation rate (see Appendix for a more rigorous argument) 



F = 37H0 + 3I3H, 



(11) 



where the parameter 7 (like /3) lies on the interval [0,1]. 
As we shall see, this scenario is compatible with the basic 
observations listed in the introduction even for P — 0. 
Inserting Eq. pT|) into ([8]) one finds 



H+-H^ 1-/3 



H 



0, 



whose solution reads 
Hit) = Ho 



1-/3; ( 
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and by integrating the above expression we obtain a big- 
bang solution for the scale factor 
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where Rq and Hq are the present day values of R{t) and 
H{t), respectively. In the limit 7^0, the above expres- 
sion reduces to 



R{t) = i?o 



2(l-/?)^ot 



(15) 



which is the model discussed in Refs. [26|, |27[, and as 
should be expected the Einstein-de Sitter cosmology is 
recovered for /3 = 0. 
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FIG. 2: The age of the Universe as a function of the 7 param- 
eter and some fixed values of /3. Note that ages great enough 
are obtained even for /3 = (bottom sohd line). For a given 
value of 7, the effect of the /9 parameter is to increase the age 
of the Universe. 



exactly the same value predicted by the 'cosmic concor- 
dance' (ACDM) model from WMAP3 and complemen- 
tary observations [3|. In the limit 7 — > one obtains 
HqIq = 2/3 as should be expected. The influence of the 
/3 parameter is apparent from Figure 2, namely, it in- 
creases the age of the Universe for a given value of 7. 

At this point, it is interesting to discuss in what sense 
this simple CDM scenario with creation behaves like an 
irreversible process. Adiabatic matter creation means 
that the total entropy S increases, but, the specific en- 
tropy (per particle), a = S/N, where N is the corre- 
sponding number of particles, remains constant |23.l23|. 
Quantitatively, a = implies that 
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Hence, due to the creation processes {N > 0), the uni- 
verse does not expand adiabatically as happens in the 
standard CDM model. Besides, since up to a constant 
factor one has N = nR^, by inserting Eq. (fTTjl into ([4]) a 
straightforward integration yields 



Ro 



(19) 



Further, from Eq. ([TH]), S" = So{N/No), and using the 
above expression one may write the entropy of the CDM 
particles like 



In Figure 1 we display the behavior of the scale factor 
as a function of time. All the models start their evolution 
from the initial singularity (i?(0) = 0). It is worth notic- 
ing that the 7 parameter does not contribute at early 
times. Actually, for H ^ Hq only the /3 parameter ap- 
pears in the equation of motion p2|) . 

Now, by taking H = H^ in Eq. ^^ or i? = i?o in 
P^ , the following expression for the age of the Universe 
is readily obtained 
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(16) 



which for 7 = reduces to Hoto ~ 2/3(1-/3) as expected 
(see [H,!!!]). 

In Figure 2 we show the age parameter as a function 
of 7 and some particular values of (3. The solid black line 
yields the age of the Universe as a function of 7 when /3 
is zero. In this case, 
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Note that ages great enough are obtained even for /3 = 0. 
In particular, for 7 ~ 0.6 the age parameter is -ffo^o ~ 1; 
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(20) 



where So is the present entropy of the CDM fluid. Note 
that if 7 = /3 = the standard conserved quantities are 
recovered. 



IV. DECELERATING PARAMETER, 

TRANSITION REDSHIFT AND SUPERNOVA 

BOUNDS 

To begin with, we first observe that by combining Eqs. 
d?]) and pT|) . the decelerating parameter reads 



^=2 



l-3/3-37f 



(21) 



so that for 7 = the value of q remains constant as 
remarked earlier. Now, by eliminating the time from Eqs. 
dni) and m]), and using that R = Ro{l + z)~\ one 
obtains the Hubble parameter in terms of the redshift 



H{z) = Ho 
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(22) 
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FIG. 3: a) Effect of the /3 parameter on q(z). For all curves 7 was taken to be zero (see Eq. (I23|l '). Note that /3 = 1/3 is a 
critical value for which g = 0. For /3 < 1/3 the possible values of q(z) are always constant and positives while for /3 > 1/3 
they remain constant and negatives in the course of the expansion. There is no transition redshift in this case, b) Effect of 
the 7 parameter on q(z). In this case, the /3 parameter has been taken to be zero. The creation of CDM particles is negligible 
at high redshifts. Due to the particle creation at redshifts of the order of a few occurs a transition from a decelerating to an 
accelerating regime. 



and inserting this result into (|2T|) it follows that 



.(^) = \ 



(1 - 3/3)(l ~ 7 ^ /3)(1 + 2)l(i-'3) - 27 
(l-7-/?)(l + z)i(i-'3)^^ 

'(23) 
(1 - 3/3)/2, while 



For 7 = 0, this expression yields q 
for /3 = we find 
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In Figure 3 we display the decelerating parameter as a 
function of the redshift as given by the above expres- 
sions. As remarked earlier, the existence of a transition 
redshift at late times depends exclusively on the 7 pa- 
rameter (compare Figs. 3a and 3b). 

A simple relation uniting 7, /3 and zt can be determined 
by taking g = 0. As one may check, Eq. (|23p implies that 
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For /3 = the above expression reduces to 

(1 + ^0^ 
^ 2 + (1 + ^0^' 
and the age of the Universe can be rewritten in terms of 
the transition redshift. One finds, 
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A. Constraints from SNe la Observations 

Let us now discuss the constraints from distant type 
la SNe data on the class of CDM accelerating cosmolo- 
gies proposed here. Since Hq can be determined from 
the Hubble Law and ^Im = 1, the model has only two 
independent parameters, namely, 7 and f3 (see Eq. ([^^ 
for H{z)). 

The predicted distance modulus for a supernova at red- 
shift 2, given a set of parameters s, is 



IJ.p{z\i 



M = 51ogdL + 25, 



(29) 



where m and M are, respectively, the apparent and ab- 
solute magnitudes, the complete set of parameters is 
s = (i/o,7,/3), and d^ stands for the luminosity distance 
(in units of megaparsecs) , 
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with x' — -4-^ ^ (l + z) ^ being a convenient integration 



variable, and H{x; s) the expression given by Eq. (|22p . 

In order to constrain the free parameters of the model 
consider now the latest sample containing 182 Super- 
novas as published by Riess and coworkers [2l- The best 
fit to the set of parameters s can be estimated by using 
a x^ statistics with 
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where ^p(z|s) is given by Eq. ^^, fil{z) is the extinc- 
tion corrected distance modulus for a given SNe la at 
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FIG. 4: a) Residual magnitudes (relative to an empty model) of 182 Supernovae data from Riess et al. (2007) and the 
predictions of the accelerating CDM model for several values of the pair (7,/3). For comparison we display the Einstein-de 
Sitter (bottom curve) and the ACDM (solid black curve). Note that the curve for 7 = 0.65 and /3 = (only one free parameter) 
is very close to the one of cosmic concordance (ACDM) scenario, b) The 7-/3 plane for a flat CDM model with gravitational 
particle creation obtained from the same sample. It should be stressed that the Supernova data can be fitted with just one free 
parameter (7) which is responsible for the transition at late times (see Fig. 3b and comments on the main text and Appendix 
B). 



Zi, and Ui is the uncertainty in the individual distance 
moduli. By marginaUzing on the nuissance parameter h 
{Ho = lOOhKm.s-^.Mpc-^) we find 0.21 < 7 < 0.75 
and < /3 < 0.46 at 95% of confidence level. The best fit 
adjustment occurs for values of 7 = 0.7 and (3 = with 
Xmin — 175.8 and i^ = 180 degrees of freedom. The re- 
duced x^ = 0-98 where (xf = Xmin/'^)^ thereby showing 
that the model provides a very good fit to these data. 



V. CONCLUSION 

In this paper we have proposed a flat cold dark matter 
cosmology whose late time acceleration is powered by an 
irreversible creation of CDM particles. In our scenario 
there is no dark energy, and, as such, the so-called co- 
incidence problem is also absent. It should be stressed 
that Hq does not need to be small in order to solve the 
age problem. Further, the transition from a decelerating 
to an accelerating regime at late times happens even if 
the matter creation is negligible during the radiation and 
considerable part of the matter dominated phase (this is 
equivalent to take /? = in all the expressions). There- 
fore, like in flat ACDM scenarios, there is just one free 
parameter, and the resulting model provides an excellent 
fit to the observed dimming of distant type la supernovae 
(see Figs. 4a and 4b). Note also that the flat model 
(J7m = 1) with creation of CDM particles proposed here 
can easily be extended to include negative (fi™ < 1) and 
positive (r^m > 1) spatial curvatures. The same happens 
with the inclusion of a small (conserved) baryonic com- 
ponent whose density parameter today is severely con- 
strained by the primordial nucleosynthesis and WMAP 
results. In this case, the value of the transition redshift 



as derived in section IV will be slightly modified. 

On other hand, the existence of such a model also 
means that the accelerating expansion does not repre- 
sent a direct evidence for a non-zero cosmological con- 
stant or, more generally, to the existence of dark energy 
as usually assumed by many authors. Naturally, new 
constraints on the relevant parameters (7 and (3) from 
complementary observations need to be investigated in 
order to see whether the matter creation model proposed 
here provides a realistic description of the observed Uni- 
verse. New bounds on these parameters coming from the 
background and perturbed equations in the presence of 
a conserved baryonic component will be discussed in a 
forthcoming communication. 



APPENDIX A: PARTICLE CREATION AND 
IRREVERSIBILITY 

In this appendix we describe how the creation pressure 
given by Eq. (O can be deduced by using the relativistic 
non-equilibrium thermodynamics. The idea is to show in 
a simplified way how an irreversible mechanism of quan- 
tum origin can be incorporated in the classical Einstein 
field equations. 

A relativistic self-gravitating simple fluid endowed only 
with gravitational matter creation is characterized by an 
energy momentum tensor T"-" , a particle current iV", 
and an entropy current S". In the homogeneous and 
isotropic case, these quantities satisfy the following rela- 
tions: 
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where (;) means covariant derivative, Pc is the creation 
pressure, n is the particle number density, F is the par- 
ticle creation rate (from quantum gravitational origin) a 
is the specific entropy (per particle) , and r is the entropy 
source. In what follows it is assumed that the particles 
spring up into space-time in such a way that they turn 
out to be in thermal equilibrium with the already exist- 
ing ones. The entropy production is then due only to the 
scalar process of matter creation (bulk viscosity has been 
neglected). Naturally, for F = we shall expect that the 
creation pressure vanishes and so also the entropy pro- 
duction. 

In the FRW background, conditions (|Aip and (|A2p can 
be written as (a dot means comoving time derivative) 



p + 3H{p + p + pc) ^0, n + 3Hn = nT. 



(A4) 



The basic aim here is to show how the second law of 
thermodynamics constrains the dependence of pc on F 
and other quantities specifying the fluid. Following stan- 
dard lines, the quantities p, p, n and a are related to the 
temperature T by the Gibbs law 



nTda = dp 



P 



dn, 



(A5) 



while the chemical potential is defined by the Euler's re- 
lation 



M 



-Ta. 



(A6) 



order to show that we consider the evolution equation 
(see section III) 



rr + Ui-^]r^ = o, 



(Bl) 



which means that the decelerating parameter {q 
—RR/B?) can be written as: 



1 
^=2 



-§ 



(B2) 



The above expression was first obtained by Zimdahl et 
al. [SQI using a different notation (see their Eq. (53)). 

Now, by taking q{zt) — in the above expression one 
finds that F = H{zt), the value of the Hubble parameter 
at the instant of transition. At low redshifts it is natural 
to take it proportional to Ho, say, F = SjHq, where the 
factor 3 is introduced for mathematical convenience and 
the constant 7 parameter, in general, depends on the 
transition redshift (see Eq. ([27]) ). Note also that the /3 
contribution can be thought as the first order correction 
of this quantity in powers of H/Hq 



= 1 



3-fHo 
and, therefore, we may write 



13 H 
7 Ho 



(B3) 



^^=2 



1 - 3/3 - 37f 



(B4) 



Now, by using equations (jA3P - 
that the source of entropy reads 



6]) it is easy to show 



r = 72<tT + Tier = — 



3Hpc pnT 



T 



T 



>0, 



(A7) 



Finally, the case of adiabatic gravitational matter cre- 
ation means that the entropy increases but the specific 
entropy a remains constant (ct = 0). Therefore, the 
above equation implies that r — nuT > with the cre- 
ation pressure assuming the form adopted in the present 
work (cf. Eq. ©) 



Pc 



P + P 
' 3H 



F. 



(A8) 



As should be expected, for F = 0, the creation pressure 
and entropy source vanish thereby recovering the perfect 
fiuid description. 

APPENDIX B: MATTER CREATION RATE AND 
THE TRANSITION REDSHIFT 

In this appendix we show a curious result, namely: the 
existence of a transition redshift, zt, at late time deter- 
mines the simplest form of the matter creation rate. In 



which is the same expression appearing in section IV (see 
Eq. dH])). For 7 = 0, the resulting scenario was proposed 
by Lima, Germano and Abramo [2^ (see also Refs. |27| ) 
while for /3 = 0, it was first discussed by Zimdahl et al. 
^30] . Clearly, the scenario proposed here is a combination 
of both approaches. Note also that only in the enlarged 
form, it may represent a possible solution to the old (and 
modified versions) of the coincidence problem (see the 
available space parameter in the (7, /3) plane as shown in 
Fig. 4b). 
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